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A variational problem for three-dimensional (3D) classical lattice models is considered with 
trial state given by two-dimensional (2D) uniform product of local variational weights. This 
approach, the tensor product variational approach (TPVA), has been applied to 3D classical 
models (the Ising and the Potts models). We consider a way of estimating the magnetic critical 
exponent (3 for the simple 3D Ising model assuming a functional form of the spontaneous 
magnetization in the off critical region, where the TPVA provides reliable data. 
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1 Introduction 

Extension of Density Matrix Renormalization Group (DMRG) f^|2||3j|4l to higher dimensions 
is still of main interest because a very slow decay in the density matrix eigenvalue spectrum 
prevents the creation of efficient renormalization group transformation |5 1. 

The numerical efficiency of the DMRG for the low-dimensional systems can be explained 
from the variational background 1 6 7 8 1 . In DMRG, the variational state is constructed by the 
product of orthogonal matrices. This type of the variational state can be defined in any dimension. 
In the three-dimensional case, Nishino et al. proposed the tensor product variational approach 
(TPVA) |9 1. It is possible to regard the TPVA as an extension of the DMRG to 3D classical (or 
2D quantum) systems. 

The TPVA has been mostly applied to 3D classical spin models such as the Ising and the 
Potts models ll9l ll0l[TTlll2l . where we considered the spatially uniform spin systems. We have 
also applied it to the 2D quantum Heisenberg model 1 1 3 1 . The TPVA is of use for the spin models 
which exhibit ordered magnetic structures such as in the 3D Axial-Next-Nearest-Neighbor Ising 
model, where the trial state is constructed from the position dependent tensors 1141 . Such a 
position dependence is also considered by Verstraete et al. for 2D quantum systems 1 15*16 17 1. 

For the models which exhibit the second-order phase transition, the TPVA tends to overesti- 
mate the critical temperature. For example, the error in Tc is about 1 % or less for the 3D Ising 
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model compared with the Monte Carlo (MC) result fl8V This is partially because the correlation 
length diverges at the criticality and the TPVA cannot treat this divergence properly due to a fi- 
nite number of variational parameters. As a result, the thermodynamic functions obtained by the 
TPVA show mean-field-like behavior around Tc. Thus, in order to obtain critical indices from 
the numerical data by the TPVA, one has to calculate the thermodynamic functions in off critical 
region, where the TPVA is more accurate. 

The main purpose of this paper is to obtain the magnetic critical exponent (3 for the 3D lattice 
models using the spontaneous magnetization calculated in the off critical region. We assume a 
scaling form for the numerical fitting of /? with respected to the calculated data. 

The paper is organized as follows: in Section 2, we define the model and discuss the vari- 
ational background of the numerical method. We propose four types of tensor product states 
(TPS) which are used to calculate the variational state. By increasing the total number of varia- 
tional parameters in the TPS, we obtain more precise numerical results. We calculate the critical 
temperature and the magnetic critical exponent for these four types of TPS in Section 3. A brief 
summary of the results obtained is summarized in Section 4. 

2 Model and variational background 

We consider the 3D Ising model on a simple cubic lattice of infinite size along all the x, y, and z 
axis coordinates. The model studied is described by the spin Hamiltonian 

= — J ^ <^i,j,k i(^i+l,j.k + (^i,j + l,k + O-'iJM+l) (1) 

with the nearest-neighbor spins cr = ±1 interacting via coupling constant J. The transfer matrix 
T connects two adjacent spin layers [ak] and [o-fc+i]. It can be expressed by product of the 
Boltzmann weights 

T[ak\ak+i] = llW^^{ak\ak+i} , (2) 

id 

where we abbreviated the notation gathering four spins into one 

Wk} = {'^id,k <^i+ld,k Cri,j + l,k (^i+l.j + l.k) ■ (3) 

We have employed the local Boltzmann weight of the interaction round-a-face (IRF) type which 
is defined as 



w!^]WkWk+i} = exp 



J , 



+ 0'id + l,kO'iJ + l,k+l + 0-i+l,j + l,kO'i+l,j + l,k+l + M<7 ,k + 0'i+l,j.,kC^i+l.j + l,k 
+ cri+l,j + l,A;fij + i,fe + (Ji^j + l^k<yi.j,k + 0'j,j,fe+lO'j+l j,fc+l + <yi+l,j,k+l<yi+l,j + l,k+l 



+ 0'i+l,j + l,A;+lf + + '^i.j + l.k+l'^idM+i) 



(4) 



where fee is the Boltzmann constant and T is temperature. Due to simphcity, we set J/k-Q = 1 
in what follows on. 
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Fig. 1 . Graphical representation of the four-spin local variational weight l^j (left) which is used to construct 
the trial state 'i/ (right). Notice that we have dropped the subscript k from the spins for simplicity. 



ai+2,j+l 




Fig. 2. The six-spin local variational weight Vij (left) and the trial state ^ (right). The weights Vij overlap 
one another only in the x axis direction (denoted by the subscript i). 



The variational partition function per layer for the given transfer matrix T has the form 

Purpose of the TPVA is to maximize Eq. Q by a proper approximation of the trial state ^. 
We assume that 4* can be written by the product of identical local variational weights V. This 
approximation is often referred to as the tensor product state Ill2lll3in"5lll6lll7l . 

In the following, we consider four candidates of ^I^, where the difference is in the definition 
of the local weights V and their connection. 

(i) In the lowest approximation, the single local weight V consists of four spins. Hence, we 
have 2"^ = 16 adjustable parameters and V graphically represents a square-shaped object 
as shown in Fig. [2 Taking product over all the uniformly given local weights V, the trial 
state yields 

^'[O'fc] = Cri+lj-,fc, CTij + i^feCTi+ij + i^fc). (6) 



(ii) A way to increase precision of the variational calculations is to enlarge the size of the local 
weight V up to six spins, see Fig. |2] Then, we let the weights V partially overlap one 
another defining the trial state as follows 

^'[O'fc] — Y\_ ^i,j{'^i,j.k,<^i+l,j,k,<7i+2J,kO'i,j+l,k,Cri+lJ + l^k,<^i+2,j+l,k)- (7) 
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Fig. 3. The local variational weight Vi,j of the IRF type with the auxiliary variables ^ (left) and the trial 
state ^ (right). The dark ovals represents the auxiliary variables to be summed up. 




Fig. 4. The local variational weight Vi^j of the vertex type with the auxiliary variables ^ (left) and the trial 
state ^ (right). The dark cubes are the summed up auxiliary variables. 



The number of the adjustable parameters is 2^ = 64. 

(ill) Further increase of the free parameters leads to a better approximation of the trial state (for 
details, see Ref. (TT\). Therefore, we introduce the so-called auxiliary variables which are 
included in the local weight V 

* K] - EH^ f?"'' 7'^"^' -^+1,.+!,'^-^ , (8) 

where ^ denotes the auxiliary variable which can be in one of n states (1, 2, n). The sum 
runs over all the auxiliary variables as depicted in Fig.|3l Thus, there are {2nY adjustable 
parameters in total. By setting n = 1, the case (i) is satisfied. 

(iv) Finally, we consider the vertex-type representation 1 19 1 of the 3D Ising model where the 
auxiliary variables ^ can be naturally included. The trial state is then given by 

^ ij V4ij,fc 4i+ij,fe 4i+i,j+i,/c/ 



and is graphically represented in Fig.|4] There are 2n'^ parameters in total. 
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Fig. 5. Left: dependence of the spontaneous magnetization (cr) on temperature T for selected variables n 
and m as listed in Table 1. Right: the plot of the effective critical exponent /3eff versus relative temperature 
t for the same variables as shown on the left graph. The total number of the adjustable parameters is shown 
in the parenthesis. 



We treat the local variational weights V with different number of adjustable parameters (up 
to 256) by the TPVA. In fact, not all of them are independent due to the symmetries of the 
model. The main advantage of the numerical algorithm TPVA is that we have derived a self- 
consistent equation (based on the a generalized eigenvalue problem) which iteratively tunes all 
the adjustable parameters regardless of initial choice. Numerical details of the self-consistent 
improvement for the trial state with explanation of the TPVA are reviewed in Refs. l^ llOlfTT] 

Having written the trial state "if and the transfer matrix T in the product forms, we can 
accurately calculate the numerator of Eq. ^ 

(*|r|vi')= l[^^A'^k}w^JWkWk+l}v,,,Wk+l} (10) 

[ak],[crk+l] i,i 

and also the denominator 

using the renormaUzation techniques. In particular, we use both the DMRG and its variant the 
Corner Transfer Matrix Renormalization Group | 4|ll9ll20l . 

3 Results 

We calculate the spontaneous magnetization per spin site using the relation 
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type of TPS 


lattice type 


' C 


'-'1 ' L ''^ J 






1 


16 


(i) 


IRF 


4.570 


1.30 


0.357 


9.58 


1 


10 


(ii) 


IRF 


4.554 


0.94 


0.348 


6.81 


2 


5 


(iii) 


IRF 


4.550 


0.85 


0.338 


3.74 


2 


16 


(iv) 


vertex 


4.533 


0.47 


0.332 


1.90 


3 


12 


(iv) 


vertex 


4.525 


0.30 


0.327 


0.37 



Tab. 1. List of the critical temperatures Tc and the magnetic critical exponents /? computed for the 3D Ising 
model for various kinds of approximation. The relative errors fx and are computed with respect to the 
results of the MC simulations 6181. 



for all the four kinds of approximation discussed above. The corresponding graph is shown in 
Fig-Elon the left. The number of the states kept in the DMRG's block-spin variable is denoted 
by the letter m. Details which concern the meaning of this m-state variable can be found in 
Ref. yj|2||3||4l. The Table 1 summarizes the critical temperatures Tc obtained by the TPVA for 
the several selected variables n and to. The MC result (Tmc — 4.5116) is used as a reference 
to evaluate the corresponding relative error £t = (7c — Imc)/7mc- Precision of our results 
improves with the increasing number of the adjustable parameters in V (for both the IRF and the 
vertex type lattice models). 

We calculate the magnetic critical exponent /? which describes the vanishing of the sponta- 
neous magnetization 

^^^^(^)'"*' ^^^^ 

when approaching the critical temperature from the feiTomagnetic phase. We analyze the 
effective critical exponent /3off as a function of the relative temperature t I21II22I 

/3eff(t) = TT^l0g(cT(t)). (14) 

dlog(t) 

Figure |5l (on the right) displays dependence of the effective exponent /3off on the reduced 
temperature t. It is evident that for t < 0.1 the linear dependence of the exponent /3off is broken 
and /Jeff tends to the mean-field value equal to 1/2. 

In order to calculate /?, we assume dependence of the spontaneous magnetization (a) on the 
relative temperature t in the series 

(f7(i)) - 7i'3'=«(*) - ^t(/3+«it+"2t^+«3t^ + -) (15) 

with unknown constants 7, /3, and on for i = 1, 2, . . .. Neglecting the terms of the second-order 
and higher in the effective exponent /3off, we have a linear dependence. The linear dependence 
is well satisfied in the case of the 2D Ising model 1221 . For simplicity, we also assume the same 
linear dependence for the 3D Ising model 1211 . in particular, 

(a(i')) = 7^' ^'"^"-'^ = It' ^+°"'' (16) 
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C> n=^, m=16, P=0.357 IRF(16) 

0.5 : A n=1. '"^lO. P=0.348 IRF(64) 

'• V «=2. "=05, p=0.338 IRF(256) 

0-4 □ «=2, m=16, p=0.332 vertex(32) i 




0.1 r 



0.0 ^ ^ 

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 

f = {T-1)/Ta 



Fig. 6. The extrapolation of tlie effective critical exponent /3eff (f') down to t' =0 using Eq. (16}. The 
plotted numerical data are not affected by the mean-field effect. The lowest full line is obtained from the fit 
of the MC simulation ['23'l. 



with t' = (Tc - T)/Tc where is obtained by the TPVA as listed in Table 1. Note, that 
Pca{t' = 0) = 1/2 is the consequence of the mean-field behavior. Therefore, we use the linear 
extrapolation of the effective critical exponent /3cff (i') = /3 + ctit' for t' ^ Q within the off- 
critical region in order to find the constants and ai. We identify the constant /3 with the 
critical exponent after taking the limit t' from the linear fit. The results are depicted in 
Fig- El Justification to carry out the linear extrapolation is supported by the plotting of the MC 
result 1 23 1 in the same graph which gives the nearly linear behavior The critical exponents /3 
obtained for various variables n and m are listed in Table 1 and compared with the MC result 
/3mc = 0.3258(14) [241. 

The error in determining (3 from the extrapolation procedure is deduced from the following 
consideration. Assume a small deviation e from the critical temperature Tc obtained by the 
TPVA, i. .e, = Tc + £. Table 2 shows that for a given e, the relative error in Tc affects the 
relative error in (3 of about ten times. 



4 Conclusion 

We have applied the TPVA to the simple 3D Ising model and calculated the critical temperature 
from the spontaneous magnetization. Using the scaling form, we obtained the magnetic critical 
exponent /3 from the spontaneous magnetization in the off critical region. The way of estimation 
of the critical index is of use for future applications of the TPVA to various 3D classical lattice 
models. 
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T'(e) 




3 




-0.010 


4.515 


0.22 


0.320 


2.1 


-0.005 


4.520 


0.11 


0.323 


1.1 


0.000 


4.525 




0.327 




H-0.005 


4.530 


0.11 


0.330 


1.1 


H-O.OlO 


4.535 


0.22 


0.334 


2.2 



Tab. 2. The relative errors and £^/3(e) calculated for given deviations e from the T^. This example is 
demonstrated on the data with the variables n = 3 and m = 12. 
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